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CN ■ Abstract 

j^ , Context. Theoretical studies of collapsing clouds have found that even a relatively weak magnetic field may prevent 

k^ ■ the formation of disks and their fragmentation. However, most previous studies have been limited to cases where the 

^—{ ' magnetic field and the rotation axis of the cloud are aligned. 

Aims. We study the transport of angular momentum, and its effects on disk formation, for non-aligned initial configu- 
rations and a range magnetic intensities. 

Methods. We perform three-dimensional, adaptive mesh, numerical simulations of magnetically supercritical collapsing 
^' dense cores using the magneto-hydrodynamic code Ramses. We compute the contributions of all the relevant processes 

f^ I transporting angular momentum, in both the envelope and the region of the disk. We clearly define what could be 

L^ , defined as centrifugally supported disks and thoroughly study their properties. 

. . Results. At variance with earlier analyses, we show that the transport of angular momentum acts less efficiently in 

f~| ' collapsing cores with non-aligned rotation and magnetic field. Analytically, this result can be understood by taking into 

O ,' account the bending of field lines occurring during the gravitational collapse. For the transport of angular momentum, 

I , we conclude that magnetic braking in the mean direction of the magnetic field tends to dominate over both the gravi- 

O ■ tational and outflow transport of angular momentum. We flnd that massive disks, containing at least 10% of the initial 

.y ' core mass, can form during the earliest stages of star formation even for mass-to-flux ratios as small as three to flve 

^2 ' times the critical value. At higher field intensities, the early formation of massive disks is prevented. 

, ^, , Conclusions. Given the ubiquity of Class I disks, and because the early formation of massive disks can take place at 

moderate magnetic intensities, we speculate that for stronger fields, disks will form later, when most of the envelope will 

p^ ' have been accreted. In addition, we speculate that some observed early massive disks may actually be outflow cavities, 

►^ ' mistaken for disks by projection effects. 

v^ ' Key words. Magnetohydrodynamics (MHD) - Stars: low mass, formation 

o\ ■ 

• r 1. Introduction as stressed by these authors, these estimat es depend on the 

^Lj assumptions made regarding the env elope (|j0rgensen et al.l 

OQ ' The formation of protostellar disks plays a central role in ' "' ™' " ™' ''' 

^H ; the context of star and planet formation. Protostars prob- From a theoretical point of view, it has been shown that, 

'^ ■ ably g row by accret ing material from protostellar accretion in the absence of a magnetic field, disks grow from small 

•!-( ■ disks (|Larson|[2003l ). and these dis ks, at lat er stages, are radii and masses by angular momentum conservatio n dur- 

X ■ the natural progenitors of planets (iLissaued 11993') . While i^S the coUapse of prestellar dense cores (e.g. Terebey_eialJ 

J^ ' observations of circumstellar disks around late young stel- |1984| ) . 

lar objects (YSO), from Class I to T Tauri stars, are well- However, observations infer that cores are mag- 
established (jWatson et alll2007[ ). it is still unclear when cir- netiz ed and typically slightly super-critical (jCrutcheil 
cumstellar disks form during the early collapse of prestellar Il999l ). that is to say the mass-to-flux ratio, M/$, is 
dense cores and the early Class phase, and what their ini- comparable to a few times larger than its critical value 
tial properties are (mass, radius, magnetic flux, and tem- ~ l/(27rG^/^). Theoretically, the presence of a magnetic 
perature). For these embedded sources, direct observations field of such intensities has been found to modif y sub- 
are indeed more difficult than for YSOs, since disk emis- stanti ally the collapse JAll en et al. 2003; Machid a et al.l 
sion is difficult to dis tinguish from the envelope signature 120051 : iFromang et al.l l2006l i and in particular the for- 
(jBelloche et al.ll2002fl. even with a relatively high spatial mation of disks. Multidimensional simulations using 



cd 



resolution (50 AU. lMaurv et al.l20TOI ). However, other stud- different numerical techniques {e.g. grid-based in 2D or 



ies observing at lower resolution (about 250 AU) and with- 3D (including adaptive mesh refinement, AMR), smooth 

out resolving the disks, infer from detailed emission model- particle hydrodynamics, SPH, codes) have shown that the 

ing the presence of disks as massive as one solar mass, cor- efficient transport of angular momentum through magnetic 

responding to about 12% of the envelope mass. Although braking may suppress the formation of a centrifugally 
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supported disk, e ven at relatively low magnetic intensities 
(li < 5 - 10, 'Mellon fc Li| 120081 : iPrice fc Bat^ 120071: 
iHennebelle fc Fromang. 2008] . Similar conclusions were 
reached by I G alii et al] (|2006f ). who performed analytical 
studies of magnetized collapsing cores. 

To illustrate the problem, we can estimate the 
strength of the magnetic field, parametrized by /i = 
(M/$)/(l/(27rG-'^/^)), for which efficient magnetic braking 
is expected. We consider the simple magnetic braking time, 
Tbr, defined as that taken by a torsional Alfven wave to re- 
distribute angular momentum from the inner to the outer 
regions of a cloud. This can be most naturally expressed as 



Tbr 






(1) 



where Zd is a characteristic scale-height, va is the Alfven 
speed, Va = ^zl \f^^^ P is the characteristic density, and 
Bz the vertical component of the magnetic field. This mag- 
netic braking time should be compared to the characteristic 
rotation time of the central region of the cloud, where a disk 
would form if the braking were not strong enough. This can 
be written as 

Trot ~ , (2) 

where rd is the disk radius, Vip the Keplerian rotational 

velocity, f^ = ^jGMjrd-, and the mass M = 2nr^ZdP- 
The ratio of these two timescales then gives 



Tbr 
Trot 



z, 



1/2 



zp, 



1 ( ^d\^ Moff 



(3) 



Where n^s is an "effective" /z for the dynamically collapsing 
structure. In general, /Zcff is smaller than the initial /x be- 
cause only a fraction of the mass contracts along the field 
line and should be used in the estimate. We recall that in a 
disk Zd < rd] the ratio of times can then be approximated 
as 



Tbr 
Trot 



< 



To"' 



(4) 



This estimate shows that for fi^s ^10 magnetic braking 
should be sufficiently efficient to remove a significant frac- 
tion of angular momentum from the inner region of the 
cloud, and thus greatly affect disk formation there. 

Recent studies have attempted to avoid this magnetic 
braking catastrophe by invoking non-ideal magnetohydro- 
dynamics (MHD) to effectively remove magnetic flux from 
the collapsing core. Ambipolar diffusion appears to be in- 
efficient: by diffusing the magnetic field out of the central 
regions, it allows instead the build up of a strong mag- 
netic field over a small circumstellar region. In this so-called 
ambipolar diffusion-induced accretion shock (Mellon & Li' 
[2009; Li et al. 2011), the magnetic braking is greatly en- 
hanced and can efficiently prevent the formation of rota- 
tionally supporte d disks. The effects of Ohm ic dissipation 
remain uncertain. iKrasnopolskv et al.l ()2010[ ) claim that an 
enhanced resistivity of about two to three orders of mag- 
nitude larger than the classical value is required although 
the physical origin of such a resisti vity remains unclear . 
However, with a classical resistivity, iDapp fc Basul (J2010h 



and lMachida et aTl (|201lf ) find tiny disks of radius the order 
of a few tens of AU, which can grow larger at later times. 
The discrepancies between these different studies of Ohmic 
dissipation might b e caused by their differen t initial con- 
ditions. In addition, ISantos-Lima et all ()2011h investigated 
the effect of turbulence: they argued that an effective tur- 
bulent diffusivity (of t he same order of ma gnitude as the 
enhanced resistivity of IKrasnopolskv et al.l 20101 is suffi- 
cient to remove the magnetic flux excess and decrease the 
magnetic braking efficiency. 

Most previous simulations have been performed in a 
somewhat idealized configuration, where the magnetic field 
and the rotation axis are ini tially align ed. As emphasized h i 
IHennebelle fc Ciardil (|2009f ) (see also iPrice fc Batell2007D . 
the results of the collapse depend critically on the ini- 
tial angle a between the magnetic field B and the rota- 
tion axis (which is the direction of the angular momentum 
J). In particular, magnetic braking has been found to be 
more efficient when the magnetic field is initially aligned 
with the rotation axis, rather than when it is not. This 
is somewhat at odds with the t heoretical conclusions of 
iMouschovias fc Paleologoul (jl979t ) , as we discuss in section 
3. 

Fo llowing th e previous studies of Hennebelle fc Ciardil 
(|2009f ) (see also'Ciardi fc Hen nebellell2010r ). we investigate 
in detail the transport of angular momentum, and the ef- 
fects of magnetic braking in collapsing prestellar cores with 
aligned and misaligned configurations (a between 0° and 
90°). 

The plan of the paper is as follows. Section 2 provides 
a general description of a collapsing core. Analytical re- 
sults describing magnetic braking are discussed in section 

3, where we show that in a collapsing core where the field 
lines are strongly bent, magnetic braking is more efficient 
when the rotation axis is parallel to the direction of the 
magnetic field, than when it is perpendicular. In section 

4, we present our numerical setup and initial conditions. 
The numerical results are presented in the section 5 and 6, 
where we focus first on the physical processes transporting 
angular momentum, and then on the physics and properties 
of the disk. Section 7 concludes the paper. 

2. Collapse, magnetized pseudo-disks, and 
centrifugally supported disks 

The gravitationally driven collapse of a magnetized core 
proceeds from an initially spherical cloud that tends to 
flatten along the magnetic fleld lines, leading to the forma- 
tion of an oblate overdensity, the pseudo -disk. Pseudo-disks 
are magnetize d, disk-like structures (jGalli fc Shu! 119931 : 
iLi fc Shulll996 'l that are not centrifugally supported. Unlike 
centrifugally supported disks, they have no characteristic 
scale and are in a sense self-similar. In this paper, we em- 
phasize the role of pseudo-disks as the place within the 
prestellar core where magnetic braking takes place. Figure 
[1] shows a slice in the equatorial plane (left panel) and along 
the rotation axis (right panel) of a dense-core collapse calcu- 
lation, for jj, — 5,a ^ 45°. In the following, we loosely deflne 
a pseudo-disk as the structure with a density n > 10^ cm""^, 
as can be seen in the right panel in Fig. [T] A more precise 
definition is given later. 

After the isothermal phase of the protostellar core col- 
lapse, an adiabatic core (the first Larson's core) with a den- 
sity > 10^° cm^'^ and a radius of about 10-20 AU forms in 
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Figure 1. Slice in density in the equatorial plane and a plane aligned with the rotation axis for /i = 5,a = 45° 
t = 24290 yr. The contours show levels of density on a logarithmic scale (n > 10^, 10^, 10^, 10^ and 10"'^° cm~'^). 



at 



the center of the pseudo-disk. This is the central object in 
Fig.[TJ We do not treat the formation of the protostar itself. 

The subsequent build-up of a centrifugally supported 
disk critically depends on the transport of angular momen- 
tum in the cloud. In contrast to pseudo-disks, which again 
are only geometrical overdensities, disks are rotationally 
supported structures formed around collapsing adiabatic 
cores. Unlike pseudo-disks, centrifugally supported disks 
possess a characteristic scale, namely the centrifugal radius. 
We discuss extensively their formation and properties in the 
next few sections. In Fig. [U the disk corresponds roughly 
to the gas with a density > 10® cm""^, and a typical radius 
of 100-200 AU (on the left panel). 

At the same time, outflows are launched in the di- 
rection of the r otation axis for all angles a < 80° 
(jCiardi fc Henneb ellc 2010). In Fig.[Tl the bipolar outflows 
are shown with an extent of about 2000 AU, and also their 
associated magnetic cavity (with a density below 10^ cm~'^, 
inside the outflows). 

The adiabatic core, the disk, and the pseudo-disk are 
embedded in an envelope, with a density of between 



3. Analytical study 

Before presenting our numerical results, we develop ana- 
lytical estimates of the braking timescales for the two ex- 
treme cases of an aligned rotator, where B and J are ini- 
tially parallel, and a perpendicular rotator, where B and 
J are initially perpendicular. The main result is that, in 
a collapsing core, magnetic braking is more efficient for an 
aligned rotator, and that disks should then form more easily 
in perpendicular rotators. 

T his result contrasts somewhat with the c lassical analy- 
ses of lMouschovias fc Paleologoul (|1979l Il980l) , who showed 
that in a simple geometry with straight-parallel field lines, 
magnetic braking is more efficient for a perpendicular ro- 
tator. However, in a collapsing core, the gravitational pull 
strongly bends the magnetic field lines, which are frozen 
in the gas, toward the center of the cloud. As suggested 
bv iMouschoviasI ()1991h . the braking efficiency should then 
increase. As we show below, using more realistic assump- 
tions that are appropriate to a collapsing prestellar core, 
the braking time for a perpendicular rotator is longer than 
for an aligned one. 



Here we focus on cores, in contra st to the classical work 
of lMouschovias fc Paleologoul (|l979l) , which was applied to 
clouds, although it does not modify the analysis. 

3.1. Magnetic braking timescales 

3.1.1. Aligned rotator {a = 0) 

Before investigating the more complex case of a core with 
strongly bent field lines, we first recall the classical anal- 
ysis performed by iMouschovias fc Paleologoul (jl979t) . The 
important point here is that in their analysis the field lines 
are straight and parallel. We first consider an aligned rota- 
tor consisting of a core of mass M, density pc, radius Re, 
and half-height Z, surrounded by an external medium - the 
envelope - of density Pext- The core has an initial angular 
velocity Q, and the magnetic field B is uniform and parallel 
to the rotation axis. A magnetic braking timescale, r|| , can 
be defined as the time needed for a torsional Alfven wave 
to transfer the initial angular momentum of the core to the 
external medium 



PoxtWA,Gxt'''|| ~ PcZ. 



(5) 



Using the expression for the Alfven speed in the external 
medium, WA.oxt = B/y/inpext , together with expressions 
for the mass of the core, M ~ 2-KpcR^Z, and the magnetic 
flux through it, $b ^ ttR^B, one obtains (jMouschoviasJ 
Il985h 

'^" ^ VPcxt/ *b' 



(6) 



This demonstrates in particular that the magnetic braking 
timescale depends only on the initial conditions, namely 
the density of the external medium, pcxt, and M/^b, the 
mass-to-flux ratio of the core. 

3.1.2. Aligned rotator {a — 0) with fanning-out 

We now consider an aligned rotator that is contracting and 
whose inner part - the core - has a density pc and a ra- 
dius Re- It is embedded in a medium of density Pext- The 
magnetic field is initially uniform. However, upon contrac- 
tion the field strength increases, through a transition re- 
gion {i.e. its envelope), from its original value -Bext in the 
external medium, to the compressed value in the core, B^. 
Magnetic braking efficiently slows down the core rotation 
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Figure 2. Schematic view of a collapsing core in aligned 
configuration with field lines that fan-out. B is the magnetic 
field, Re the radius of the core, Rq the initial radius, Z the 
half-height of the transition region, pc the density of the 
core, and pext the density of the external medium. 

if this braking sets into co-rotation an amount of matter in 
the envelope with a moment of inertia equal to that of the 
core. Since the angular momentum in the envelope must 
be approximately equal to the momentum of the core, and 
assuming a co-rotation of the field lines, we find that 



Oi?o T^Ro PcxtVA,cxtTn Jo 



nRl ttRI PcZ, 



(7) 



(8) 



where Rq is the initial radius of the flattened core, Z its 
half- height, t'A,ext the Alfven speed in the external medium 
(^A,ext = -Bext/V47rpcxt), and ry the magnetic braking 
timescale for the aligned configuration. We thus obtain the 
magnetic braking time for the case of fan-out, t^\joi which 
is given by 

Pc Z f Re 

Poxt WA,cxt V R-O , 

This expression is identical to Eq. [SI apart from the coeffi- 
cient {Rc/RqY, whose origin is twofold. First, because the 
field lines fan-out, the volume of the external medium swept 
by the Alfven waves increases more rapidly than when the 
field lines are parallel. This accounts for a factor {Rc/Rq)'^. 
Second, as co-rotation of the field lines is assumed, fluid 
elements of the external medium lying along diverging fleld 
lines have higher speciflc angular momentum than if they 
were on straight fleld lines. This accounts for another factor 
{Rc/Ro) ■ 

Using again the mass and the magnetic flux of the core, 
M ~ 2'iTpcR^Z and <^b ^ 7ri?o-^cxt — T^R^Bg and the ex- 
press ion for WA ext , we can rewrite Eq. ([S]) as (jMouschoviasI 
119851) 



nijo 



VPcxt/ 



1/2 



M f R, 



*B \R 



(9) 



Therefore, when the fleld lines fan-out, the magnetic brak- 
ing timescale depends not only on poxt and M/^b, but 
also on the contraction factor of the core Rc/Ro- In a 
collapsing core Re ^ i?Oj and this geometrical factor 
can signiflcantly reduce the characteristic magnetic brak- 
ing timescale in an aligned rotator, the ratio of timescale 
being t\\/t\\jo = {Ro/Rcf- 



Although this analysis gives a more realistic estimate of 
the braking time in a collapsing prestellar core, it is still 
greatly simplified. In particular, it is assumed that the field 
lines are immediately set into co-rotation. Because of the 
collapse, as the waves propagate outwards, the radius i?o 
and the density pext vary continuously, before reaching ap- 
proximately constant values. 



3.1.3. Perpendicular rotator [a 
decreasing Alfven speed 



90°) with radially 




Figure 3. Schematic view of a collapsing core in perpen- 
dicular configuration. B is the magnetic field. Re the radius 
of the core, pc the density of the core, and poxt the density 
of the external medium. 



In the case of a perpend i cular rotator, the analysis 
of iMouschovias fc Paleologoul (|1979f ) considers the brak- 
ing timescale corresponding to the time needed for Alfven 
waves to reach i?j^, the radius for which the angular mo- 
mentum of the external medium is equal to the initial an- 
gular momentum of the core. In this case, Alfven waves 
propagate in the equatorial plane, rather than along the 
rotation axis, and sweep a cylinder of half-height Z and 
radius _Rj^, thus 



/5cxt(^i 



Re 



PcRt- 



(10) 



Assuming further that the magnetic fleld has a radial de- 
pendence, B{r) (X r^^, so that VA{r) — va{Rc) x Rc/r, the 
perpendicular rotator magnetic braking time is then given 
by 



r± 



1 Re 



Re 



dr 



VAir) 2va{Rc) 



\ 1/2 
1 + ^) -1 

Pcxt/ 



(11) 



Using the expressions for the mass, M — 2iTpeRcZ, mag- 
netic flux, $B — ATrReZB{Rc)., and Alfven speed VAiRe), 
together with the approximation pc S> Pcxt, Eq. [11] becomes 



T_L 






(12) 
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3.1.4. Perpendicular rotator [a = 90°) with constant Alfven 
speed 

The expression in Eq. (J12p was derived assuming a magnetic 



field B{r) oc 



and consequently an Alfven speed that 



decreases with radius. However, the field lines are twisted 
because of the rotation of the core and are not purely ra- 
dial. It can easily be inferred from the divergence-free con- 
straint on the magnetic field, and from the simulations, that 
the structure of the magnetic field is more complex, and is 
not radial, as shown schematically in Fig. [31 It is there- 
fore unlikely that the Alfven speed drops with radius as 
1/r, and collapse calculations indeed show that the Alfven 
speed remains roughly constant in the dense cores (e.g. 
iHennebelle eTall l201lh . 



given by 



T±., 



^ dr 



In this case, the braking time is 



R± — Re 



VA 



Re 

VA 



Pc 
Pcxt 



1/4 



- 1 



Since Pc/pcxt ^ 1, the braking time then becomes 

:Plit M 



'''^,CVA ^ 4-^7?- 



Pc 



3/4 $^ ■ 



(13) 



(14) 



This braking time is shorter than that obtained for a radi- 
ally decreasing Alfven speed (cf. Eq. [T^ , their ratio being 



T-l _ 1 / Pc \ 
T^-,cvA 2 \PcxtJ 

3.1.5. Comparison of timescales 



1/4 



(15) 



in the preceding sections, we have derived four character- 
istic magnetic-braking timescales: two for aligned rotators 
consisting of one with straight field lines, rii , and one with 
field lines that are fanning-out, t\\joj ^-nd two for perpen- 
dicular rotators, namely one for a radially decreasing Alfven 
speed, T±, and one for a constant Alfven speed, t±_cva- 
Comparing the braking timescales given in Eq. [S] and \TI\ 
gives 

1 / \ 1/2 

T± 2 VPcxt/ 



(16) 



and since pc 3> Pcxt, this leads to the conclusion that the 
magnetic braking is more efficient in the perpendicular case 
than the aligned one (Mouschovias & Palcologou 1979J3- 
This conclusion may apply to prestellar cores whose density 
is not centrally condensed, and for which the fanning-out 
is weak. It probably applies to the outer part of the cores, 
making it quite possible that during the prestellar phase, 
before the collapse occurs, magnetic field and angular mo- 
mentum may become, to some extent, aligned. Owing to 
the turbulent motions in the ISM, it is not unlikely how- 
ever that the cores have a misalignment between the rota- 
tion axis and the magnetic field. But this conclusion that 



^ Although R and Z do not appear explicitly in this expres- 
sion, we recall that the definitions of the magnetic flux are not 
the same in both cases. A more correct expression should include 
a factor <l?||/<l?x. 



the magnetic braking is more efficient in the perpendic- 
ular configuration does not apply to the internal part of 
the collapsing cores where magnetic field lines are strongly 
squeezed toward the center. 

To be more quantitative, we consider the external 
medium to be the core's envelope, and estimate the average 
Pcxt to be a few times the density of the singular isothermal 
sphere, Pcxt{Ro) oc Rq , and similarly for the core density 
Pc{Rc) fx Rc^. Using Eqs. © and (fT^ . gives for the ratio 
of timescales 

T± Rq' 



(17) 



Since Rc/Rq <C 1, the angular momentum is more effi- 
ciently transferred to the envelope in an aligned rotator 
than in a perpendicular one. This is still so when consider- 
ing the more realistic case of a perpendicular rotator with 
a constant Alfven speed. Using Eq. ^ and ([H]), the ratio 
of the timescales is 



nijo 



T-L. 



R^ 
Ro 



1/2 



(18) 



Evidently, the previous conclusion still holds in this case; 
magnetic braking is more efficient in an aligned configura- 
tion than in a perpendicular one, although the difference is 
smaller. As both configurations somehow represent extreme 
cases, we expect our simulations to have properties that are 
in-between. 

To conclude, there are four magnetic-braking timescale 
of interest, which have the following ordering tii > t± > 
T±.cvA > ''"IL/o- The first two inequalities hold in (non- 
collapsing) prestellar cores, whose density is not centrally 
condensed and fanning-out is weak. However, for conditions 
that are more appropriate to collapsing prestellar cores, it is 
the last inequality that is appropriate, and aligned rotators 
are more efficiently braked than perpendicular ones. 



4. Numerical setup and initial conditions 

4.1. Numerical setup 

We perform three-dimension al (3D) numerical simulations 
with the AMR code Ramses (|Tevssieill2002l:lFromang et all 
I2006D . Ramses can treat ideal MHD problems with self- 
gravity and cooling. The magnetic field evolves using the 
constrained transport method, preserving the nullity of 
the divergence of the magnetic field. The high resolution 
needed to investigate the problem is provided by the AMR 
scheme . Our simulations are p erformed using the HLLD 
solver (|Mivoshi fc Kusanoll2005l ). 

The calculations start with 128^ grid cells. As the col- 
lapse proceeds, new cells are introduced to ensure the Jeans 
length with at least ten cells. Altogether, we typically use 8 
AMR levels during the calculation, providing a maximum 
spatial resolution of ~ 0.5 AU. 

4.2. Initial conditions 

We consider simple initial conditions consisting of a spher- 
ical cloud of 1 Mq . The density profile of the initial cloud 

Po 



1 + {r/ray 
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Table 1. List of performed calculations 



where po is the central density and tq the initial radius 
of the spherical cl o ud is in accordance with observations 
(| Andre et al.ll2000t iBelloche et all I2OO20 . The ratio of the 
thermal to gravitational energy is about 0.25, whereas the 
ratio of the rotational to gravitational energy /3 is about 
0.03. We run 17 simulations. Various magnetization cases 
are studied: /i = 2, 3, and 5 (magnetized super-critical 
cloud, in agreement with observations, as pointed out in 
the introduction) and 17 (very super-critical cloud). The 
angle between the initial magnetic field and the initial ro- 
tation axis a is taken to be between and 90°. Table[T]lists 
all the simulation parameters. 

To avoid the formation of a singularity and mimic that 
at high density the gas becomes opaque i.e. nearly adia- 
batic, we use a barotropic equation of state 



P ' 



^2 



VPad/ 



2/3' 



where pad is the critical density over which the gas be- 
comes adiabatic; we assume that pad = 10"^'^ g cm~'^. 
When p > Pad, the adiabatic index 7 is therefore equal to 
5/3, which corresponds to an adiabatic mono-atomic gas. 
At lower density, when the gas is isothermal, P/p is con- 
stant, with Csfi ^ 0.2 km s~^. The corresponding free-fall 
time is tg ^ 12 kyr (for an initial density peak of about 
3 X 10^1^ g cm-3). 

5. Transport of angular momentum 

We analyze in detail the transport of angular momentum 
in our numerical simulations. 



5.1. Temporal evolution 

We begin by describing the temporal dependence of the 
angular momentum, in particular the specific angular mo- 
mentum, which is defined by 



^-^ 



r X pv dV, 



(19) 



where M is the mass contained within the volume V, r 
the position (with respect to the center of mass), p the 
mass density, and v the velocity. In general, we compute 
three values of J using three density thresholds corre- 
sponding roughly to the adiabatic core {n > 10^° cm 



the disk {n > 
envelope (n > 



10^ 
108 



cm ), and the densest parts of the 
cm~^). These are nested structures: 



by n > 10^ cm~^. Figure [H displays the norm of the 
specific angular momentum |J| for all considered orienta- 
tions, magnetizations (p — 17, 5, 3) and density thresholds 
(n > 10^, 10^ 10^° cm-3). Figure displays |J| for all the 
density thresholds for fj, — 5 and three different orientations 
(a = 0, 45,90°). 

Figure [5] illustrates that, as expected, the angular mo- 
mentum increases with decreasing densities, which corre- 
spond to the outer regions of the collapsing prestellar core. 

Figure [3] shows that, as matter is continuously ac- 
creted, the angular momentum increases with time, and 
it is smaller for larger magnetizations, an indication that 
magnetic braking is more efficient in transporting angular 
momentum from the inner to the outer parts of the prestel- 
lar core. 

As for the dependence on the angle a, there are several 
interesting aspects worth discussing. First of all, for the 
disk (n > 10^ cm~'^) and adiabatic core (n > 10^° cm"'^) 
the angular momentum increases with a (see Fig. [3]). 
Conversely, magnetic braking rapidly decreases with a, 
which is consistent with the prediction that the magnetic 
braking timescale for a perpendicular rotator is longer than 
for an aligned one (see section [3]). For the disk region, the 
angular momentum in the perpendicular case is indeed al- 
most three times larger than in the aligned one {p = 5). 
This proportion decreases with the magnetization: it is a 
factor of two when p = 3, and for p = 2 the values are 
comparable, albeit the angular momentum still increases 
slightly with the angle a. Therefore, in misaligned rota- 
tors and for intermediate magnetizations, more angular mo- 
mentum will be available to "build" centrifugally supported 
disks. 

Interestingly for p — 5 and 3, the angular momentum 
below a density of 10® cm~^ is independent of a, which sug- 
gests that its efficient transport only occurs in the highest 
density regions {n > 10® cm^'^). This corresponds therefore 
to a "braking region". This is no longer true for p = 2: mag- 
netic braking occurs earlier, simply because the magnetic 
field is stronger. 

5.2. General considerations 

The azimuthal component of the conservation of angular 
momentum, in cylindrical coordinates, is the starting point 
of our analysis. In conservative form, it is given by 



dt (prv^) 





r / 


Vt 


pv^v + 

V 


47r 


47rG^ 



p 






0, 



(20) 



the adiabatic core is embedded in the structure described 



where p is the density, v the velocity, P the gas pressure, B 
the magnetic field, and the gravitational acceleration g — 
V<I>, where $ is the gravitational potential. For the sake of 
completeness, its derivation is detailed in the appendix. 

The corresponding fluxes of angular momentum in this 
equation are rpv^v for the mass flow, rB^'B/Air for the 
magnetic field, and rg^g/Air for the gravitational field; they 
represent the contribution of each of those processes to the 
transport of angular momentum. In general, magnetic brak- 
ing can be efficient in both the vertical and radial directions. 
However, the outflows remove angular momentum mainly 
in the vertical direction (mass accretion occurs instead in 
the radial direction, carrying in angular momentum). The 
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Figure 4. Evolution of the specific angular momentum jj [J ^ r x pv dV) for ^u = 17 (Fig. 4(a) \, fi ~ 5 (Fig. |4(b)[ ), 

fj, = 3 (Fig. 4(c)), and /x = 2 (Fig. |4(d")] ), for three different density thresholds pcr that correspond to: n > 10^° cm~^, 
n > 10^ cm '^, and n > 10^ cm^^. 



gravitational transport of angular momentum is most effi- 
cient in the radial direction, because of the spiral arms that 
develop around the first core. 

The pressure terms (P + B^/Stt — g'^/SirG) do not 
contribute significantly to the transport of angular mo- 
mentum. In the following, to quantify the contribution of 
each of these processes, these fluxes are considered and 
compared one by one. One should obviously add all these 
terms to obtain the total angular momentum transport. 



To calculate the fltixes defined above, one needs to de- 
fine the main axis of a cylindrical frame of reference. Two 
different choices have been made below: we use either the 
inertia matrix or simply the rotation axis of the system. 

The pseudo-disk is defined as all matter with a parti- 
cle number density n > 10^ cm^'^. We then calculate the 
inertia matrix over the volume V of the pseudo-disk as 



J^ij — 



pnrj 



dV, 



(21) 
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where p is the density, and r^ the coordinates (with «,j G 
{1,2,3}, so {ri,r2,r3} = {x,y,z}) of a fluid element with 
respect to the center of mass. The eigenvectors of this ma- 
trix are the axis of the frame of the pseudo-disk, which we 
call TZp in what follows. The z-axis of the frame TZp is de- 
fined as the eigenvector associated with the minimum eigen- 
value of the matrix. It corresponds to the z-axis of the cylin- 
drical coordinates in the following analysis. As the pseudo- 
disk is essentially perpendicular to the magnetic field, the 
z-axis of this frame is close to the direction of the actual B. 

The fluxes defined above are then computed on surfaces 
at the edge of the pseudo-disk. Annuli are then defined 
to fit its surface: we consider a set of one hundred annuli 
of radius r, between and Rq, height h{r) = r/4, and 
thickness i?o/100, where Rq is about 1000 AU. The fluxes 
are computed on the cells belonging to these annuli. 

We note that we also tried to estimate the actual value 
of h{r), as a function of radius, although we found that 
the surface integrated fluxes did not change signiflcantly. 
In addition, we also considered a frame whose main axis 
is the rotation axis of the pseudo-disk instead of the main 
axis of the inertia matrix, and our conclusions remain again 
qualitatively unchanged. 

For the disk, we simply choose for the cylinder axis, the 
rotation axis of gas denser than 10^" cm"^, which corre- 
sponds to the rotation axis of the core itself. We refer later 
to this frame as TZd- We define annuli as previously, but 
restrict the analysis to a maximum radius of 400 AU since 
disks are not larger. 

The integrated vertical and radial fluxes of angular mo- 
mentum transported by the magnetic fleld are then defined 
as 



F^'^iR) 



2tc 



B4r,cP,±h{r)/2)B,{r,cP,±h{r)/2) . 

r — rdr c 

47r 



and 



-h{B.)/2 



(22) 



^ B4R,c^,z)Br{R,(f>,z) _^^_^^ 
47r 



where Bi{R,(j},z) = Bi{r — R,<l),z). We note that F^ is 
the sum of the fluxes through the faces defined by h{r)/2 
and —h{r)/2. To compare the various cases and time-steps, 
it seems appropriate to consider specific quantities. In the 



foUowing, we study F^/M and F^/M, where M is the 
mass enclosed in the volume of interest. We do the same 
for the outflows and the gravity terms. 

5.3. Transport of angular momentum in the envelope 

We flrst investigate the transport of angular momentum 
in the envelope (see section [5]), focusing on the magnetic 
braking. As we show below, this is the most efficient means 
of transporting angular momentum, in particular in the 
strongly magnetized clouds. We work here in the frame TZp. 

Figure 6(a) shows the evolution of specific radial flux 
{F^ /M) in the left panel, and specific vertical flux {F^ /M) 
in the right panel, at four different time-steps, for /i — 
5,a = 0°. 

The two panels of Fig. 



a) show that the magnetic 



braking depends on the radius: it is more efficient in the 
inner region of the envelope than in the outer region. In 
the outer part of the cloud, both components decrease with 
the radius r. In particular, _B,. and B^ drastically decrease 
outside the cavity of the outflows (see section [5]) since the 
twisting of the magnetic field lines - which generates both 
the radial and azimuthal components of the magnetic field 
- occurs essentially inside the cavity of the outflows. 



Figure 6(a) shows that the vertical component of the 
magnetic braking is larger than the radial component (by 
about one order of magnitude in the inner part of the en- 
velope) . The ratio of those two components increases with 
the radius, since the radial component of the magnetic fleld 
is almost zero at large radius. 

Figures 6(b) and 6(c) display a comparison of the ra- 
dial (left panel) and vertical (right panel) components of 
the magnetic braking for a = 45° and 90° with th e com po- 
nent s in the aligned case. The left panels of Fig. 6(b) and 
6(c) show that the magnetic braking in the radial direction 
is less efficient in the tilted cases than in the aligned case in 
the inner part of the envelope by a factor ~ 2—10 (especially 
for a = 90°). The limit of the cavity is shown by the sharp 
increase in the ratio (around ~ 200 AU for the first time- 
step, and ~ 500 AU for the later time-step). Outside the 
cavity, the ratio increases because Br decreases more dra- 
matically in the aligned case than in the misaligned cases; 
this is because the radial component vanishes initially in the 
aligned case whereas there is an initial Br in the misaligned 

cases. 

show that mag- 



The right panels of Fig. 6(b) and 



netic braking is less efficient in the vertical direction in the 
tilted cases than in the aligned case. This is true everywhere 
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Figure 6. Magnetic tra nspor t of ang ular m omentum in logarithmic scale, for ^ = 5, a = 0° (Fig. 6(a)), a = 45° (Fig. 
16 (b) I, and a ^ 90° (Fig .|6(c) [). Figure |6 (a) I represents log(i^^^/M) (left panel) and log(F„^/M) (right panel) for four time- 



steps. Figures [6(b)] and |6[cjfdisplay l(^F^ / M) / {F^^ / Mq)) and \og{{F^ /M)/{F^,^/Mo)) where i^ifoA^o and F^^oA^o 



correspond to the aligned case, at corresponding timesteps 



by a factor of ^ 2 — 5. As we see later, the vertical compo- 
nent dominates the transport of angular momentum by the 
magnetic field, it clearly shows that the magnetic braking 
is more efficient in the aligned case than in the misaligned 
cases; this is consistent with our previous analytical analy- 
sis. 



5.4. Transport of angular momentum in the disk 

In the region of the disk, magnetic braking, outflows, and 
gravitational torques all contribute to the transport of an- 
gular momentum. We here work in the frame TZd- 



5.4.1. Comparison between magnetic braking and extraction 
by the outflows 

Outflows are one of the most important tracers of star for- 
mation. From the very beginning of protostar evolution 
to the T Tauri stage, they are t hought to be launched 
by magneto-centrifugal means fiBlandford fc: Pavng Il982t 
iPudritz fc NormanI Il983l . lUchida fc Shibatal 119851) . and 
may play an important role in the efficient transport of 



angular momentum (jBacciotti et al.l [20021 ). The early for- 
mation of outflows during the collapse of dense cores 
was i nvestigated recently by 2D and 3D MHD sim ula- 
tions (jMellon fc Lill2008HHennebelle fc Fromandl2008D. and 
where t he second coll apse was included fe.g. lBaneriee et all 
2006; M achidaet an i20081. 

To study the impact of the outflows on the transport of 
angular momentum, we begin by comparing the integrated 
flux of angular momentum transported by the magnetic 
field, Finag, with that of the outflows. Font, as a function of 
time. The respective integrals are given by 



^ mag — 



An 



for the magnetic braking. 



-Fout 

for the outflows, and 



/ prv^v ■ dS 
Js 



Fin = 



/ prv^v ■ dS 
Js 



, V • dS > 



, V • dS < 



(23) 



(24) 



(25) 
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Figure 7. Evolution of angular momentum trans ported by the mag netic field and by the outflo ws within a cylinder of 
radius 300 AU and height 150 AU, for /x = 5 (Fig. [7(a)] ), /i = 3 (Fig. [7(b)] ) , and /i = 2 (Fig. [7(c)] ) . 



for the accretion flow. 

The integrals are taken over the surface S* of a cylinder, 
corresponding approximately to the disk, of radius i? ~ 
300 AU and height h ~ 150 AU, whose axis is taken to 
be that of rotation. As previously, we study the specific 
quantities Fmag/M, Fout/M, and Fin/M in the following. 

Figure [7] displays the evolution of angular momentum 
carried away by magnetic braking (left panel) , and the out- 
flows (right panel) for /i — 5,3, and 2, respectively. More 
precisely, the ratio of this quantity to the total mass en- 
closed in S is computed. The left panel of Fig. [7] shows that 
more angular momentum is carried away from the central 
region of the collapsing core for relatively small a (below 
70°) than for larger a (above 70°). This is the case for all the 
magn etiz ation s but is particularly evident for lo wer /x (Fig. 



7(b) [ and [7(c 
that the angu 



). In the right panel of Fig. 7(a) 



we can see 



ar momentum carried away by the outflows is 



comparable to that transported by the magnetic braking in 
the aligned case. When the angle a increases, the amount of 



angular momentum carried away decreases and in the per- 
pendicular case, the total angular momentum transported 
by the flow is about ten times smaller than in the aligned 
case. The suppression of the outflows with i ncrea sing a is 
clearly responsible for this decrease. Figures 7(b) and 7(c) 
show that this effect is even stronger, since the increasing 
magnetic intensity (i.e. decreasing n) reduces the strength 
of the outflows. Less and less momentum is therefore car- 
ried away by the outflows with increasing a and decreasing 
/i. 



Figure [5] displays the ratio of the flux of angular mo- 
mentum transported outward (by the magnetic field, Fmag, 
and the onflows. Font) to that transported inward (by the 
accretion flow, Fm) for /i = 5. While accretion dominates 
for a > 45°, the angular momentum both accreted and 
expelled are comparable in the other cases. This can be un- 
derstood by recalling that, in steady-state, Eq. (|^U| reduces 
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ported outward (Fmag + Fout) to the flux of angular mo- 
mentum transported inward (Fin), within a cylinder of ra- 
dius 300 AU and height 150 AU, for /i = 5. The straight 



line corresponds to {F„ 



^ out ) / -^ in — -L- 



gular momentum from the central part of the cloud (within 
a radius of 100-150 AU, where th e rati o is meaningful). 

In the misaligned cases (Fig. 9(b) and 9(c)), it is also 
clear that magnetic braking dominates in the central part 
of the cloud. 



5.4.2. Gravitational transport of angular momentum 




(a) 



to 



V- 



P^'0V - — B 



= 0, 



(26) 



where the other terms, in particular the gravitional torques, 
have been neglected (see discussion in section I5.4.2p . The 
steady-state condition therefore reduces to F^ag + Font ^ 
Fin, which is approximately the case for a < 45°. 

To get a deeper understanding of the impact of the flows 
on the transport of angular momentum, we look at the spa- 
tial distribution of the fluxes and compute the total flux of 
angular momentum transported by outflows and magnetic 
braking in concentric cylinders of constant height H. Those 
cylinders are oriented along the rotation axis of the core, 
since the outflows are approximately aligned with it. We 
consider only the mass expelled from the core {i.e. with a 
positive vertical velocity) hence only the vertical compo- 
nent of the flux, since mass is accreted mostly along the 
radial direction. Figure [5] displays the vertical flux of an- 
gular momentum transported by the magnetic field (left 
panel) and the outflows (right panel), for /i = 5 and three 
different angles (0°,45°, and 90°). The integrated flux of 
the angular momentum carrie d by the magnetic field (left 
panel of Fig. 9(a) 9(b) and 9(c) I increases with the ra- 
the outflows (its radius 
There, a local reversal of 



dius, until the limit of the cavity o: 
is from about 150 to 200 AU) 
the magnetic field usually happens, which provokes a local 
variation in the integrated flux. Outside the cavity, B^ is 
close to 0, which means that braking no longer occurs any- 
more and the integrated flux remains almost constant (the 
differential flux (5Fniag being close to zero). 

The integrated flux of angu lar momen tum carried by 
the outflows (right panel of Fig. 9(a) 9(b) and 9(c)) sim- 
ilarly increases with the radius inside the cavity. Outside 
the cavity, no outflow occurs and the integrated flux 
therefore remains a l most constant. As pointed out in 
ICiardi fc Hennebelld (J2010f ) , there is almost no outflow in 
the perpendicular case (Fig. 9(c)). Thus, almost no angu- 
lar momentum (or more precisely one o rder of magnitude 
less than in the aligned case. Fig. 9(a) I is transported by 
the outflows in this conflguration. A comparison of the two 
previous integrated fluxes conflrms that magnetic braking 
is globally more efficient than the outflows in removing an- 
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Figure 10. Gravitational transport of angular momentum 
for fj, = 17 and a = 0, 45, and 90°, at /i = 46 AU. It repre- 
sents the ratio (in logarithmic scale) of the radial compo- 
nent of gravitational (FJ/A/) to magnetic (F^/M) trans- 
ports of angular momentum. 



The integrated gravitational ffuxes of angular momen- 
tum are given by 



h/2 



-h/2 Jo 



'\^ 



R, z)g,r{r — R, z) 



AtiG 



27rM(/>dz 



(27) 

As for the other fluxes, we compute the specific quantity 
Fg/M. We focus on the radial component of these angular 
momentum transport processes because the gravitational 
transport acts mostly in the radial direction. 

Figures [10] and [11] show a comparison of the magnetic 
and gravitational braking in the region of the disk for fi = 2 
and /z = 5. It represents the evolution of the logarithm 
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Figu re 9. Magnetic tran sport -^ J rB^H/A -K ■ dS and transport by the outflows j^ J rpv^v ■ dS, ^ — 5, a = 0° (Fig 



9(a)]), a = 45° (Fig. [9(b)| ) , and a = 90° (Fig. [9(c)| ) 



of the ratio of their radial components on the surface of 
cylinders of fixed height {h — 46 AU) and a radius between 
and 400 AU. The height of the cylinder is defined to take 
into account the whole gravitational transport in the disk. 
There is no (or negligible) gravitational transport above 
and below the disk and the typical height of an hydrostatic 
disk is about 30 AU (see section [6^ : integrating over these 
cylinders gives the whole gravitational transport in the disk. 

The /i = 17 case (Fig. \W\\ corresponds to a quasi- 
hydrodynainic case, where the magnetic field strength is 
low. The magnetic braking is less efficient than in higher 
magnetization cases, and in every configuration, a centrifu- 
gally supported disk forms. Therefore, the gravitational 
transport of angular momentum is always larger (up to ten 
times larger) than the magnetic transport of angular mo- 
mentum in the disk (within a radius of 150 AU). Outside 
the disk, the gravitational transport is weaker (10 to 100 
times weaker than the magnetic transport). 

Since there is no disk in the aligned case for fJ, — 5, (Fig. 



11(a)), the radial component of the gravitational trans- 
port of angular momentum is about 10 to 30 times weaker 
than the radial component of the magnetic braking under 
150 AU and around 500 times weake r for la rg er rad ii. In 
the misaligned cases, for /i = 5 (Fig. 11(b) & 11(c) I, the 



gravitational contribution to the transport of angular mo- 
mentum gradually increases: it is 10 times weaker than the 
magnetic transport for a = 45° and of the same order of 
magnitude to 3 times weaker in the perpendicular case, 
for radius r < 100 AU. For r > 100 AU, for all the mis- 
aligned cases the gravitational contribution is between 10 
and 100 weaker than the magnetic one. On the one hand, 
the magnetic transport becomes weaker as a increases and 
on the other hand, gravity transports momentum more ef- 
ficiently in the presence of a disk, owing to density waves 
(the spiral arms) that propagate in the radial direction. 
The gravitational transport is nonetheless less efficient be- 
cause of the symmetry of the disks, which is stabilized by 
the m agnetic field, as emphasized in lHennebelle fc Tevssieij 
(|2008f ) ; less symmetric disks would transport more momen- 
tum by means of gravitational torques. The gravitational 
transport is stronger in the perpendicular case than in the 
other misaligned cases because disks are more massive. 

For other magnetizations, these conclusions hold, since 
without a disk the magnetic braking is the most efficient 
process of angular momentum transport, whereas in the 
presence of a disk, a significant - although not predomi- 
nant - fraction of the momentum can be transported by 
gravity. In the low magnetized cases {fi = 17), gravity can 
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Figure 12. Density map on a logarithmic scale in the disk 
for /i = 5,a = 90°, for three different time-steps. The disk 
grows with time as the central part becomes denser; the 
maximum radius corresponds to the edge of the spiral arms 
of the disk structure. 



transport even more momentum in the radial direction than 
the magnetic field. 



6. Disk properties 

When enough angular momentum is left in the envelope, 
a disk can form around the adiabatic core. Here, we first 
discuss in detail how to define disks, then study some of 
their properties. 



6.1. Disk formation 

We work in the frame of the disk TZd (see section lS^ . where 
the main axis of the frame is defined by the angular mo- 
mentum. Several criteria must be used to define disks. As 
we show below, a simple rotation criterion is insufficient 
to define a disk because several parts of the envelope are 
rotating but do not belong to the disk. For example, defin- 
ing the disk as all material whose rotation velocity is larger 
than a few times the infall velocity would also pick up the 
walls of the cavity of the outflows, which are rapidly rotat- 
ing. A single geometric criterion is also insufficient; disks 
are not well-approximated by cylinders. We define disks by 
employing a combination of five different criteria. As disks 
are expected to be reasonably axisymmetric, these criteria 
are defined for concentric and superposed rings in which 
density, velocity, pressure, and magnetic field are averaged. 



As disks are expected to be Keplerian, we first use a 
velocity criterion. A ring of matter should not collapse 
too rapidly along the radial direction, which implies that 
the azimuthal velocity must be larger than the radial 

velocity {v^ > /thresWr)- 

As disks are expected to be near the hydrostatic equilib- 
rium, the azimuthal velocity is larger than the vertical 
velocity {v^ > fthrosVz)- 

The central adiabatic core is also rotating but is not in 
the disk; another criterion is therefore added, to take 
into account only areas of the simulation that are ro- 
tationally supported. We thus check whether the rota- 
tional support (the rotational energy, pw?/2) is larger 
than the thermal support (the thermal pressure, Pth) 
by a factor /thrcs- 

A connectivity criterion is also used: a ring area belongs 
to the disk if it is linked to the equatorial plane. 
As discussed in section 16.41 we add a density criterion 



{n > 10 cm ^) to avoid the large spiral arms and ob- 
tain more realistic estimates of the shape of the disk. 

For the three first criteria, a value /thrcs = 2 is chosen 
below. 

Figure [12] shows the azimuthally averaged shape of the 
disk using these criteria, for /i = 5,a = 90°, and three 
different time-steps. 
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6.2. Mass 

The mass of the disk as a function of time, for different mag- 
netizations and angles, is presented in Fig. [T31 The general 
trend shows an increase in the disk mass with the angle a. 
This agrees with our previous discussions, which indicated 
that the braking time in the parallel case is shorter, lead- 
ing to a more rapid removal of angular momentum in the 
infalling envelope, thus limiting the effective mass of disks. 
At the same time, it is clear that for increasing magnetic 
field strength, thus increasing magnetic braking, disks with 
masses greater than O.OSMq are only found in misaligned 
configurations. The limiting case corresponds to a magne- 
tization of fi — 2, where the removal of angular momentum 
by the magnetic field is so efficient that the mass of rotat- 
ing gas does not exceed 0.05Mq, even in the perpendicular 
case. We note that higher resolution simulations tend to 
ind icate m asses that are even smaller than this value (see 
fig.Efb)). 

At the other end for low magnetization (/i = 17), disks 
always form with masses that increase to about O.SMq. We 
note that in the aligned case, the disk fragments, leading 
to a decrease in its mass after a time i ~ 24 kyrs. 

In the intermediate regime of magnetizations (/i = 5), 
the mass of the disk starts to increase significantly even 
for small a. For the intermediate angles (20 and 45°), disk 
masses increase within the range 0.15 — 0.2Mo. For the 
more tilted cases (70, 80°, and the perpendicular case), 
disks grow to 0.25 - O.4M0. 

For /J, = 3, disks do not form for angles < 20°. In the 
45° case and the perpendicular case, the masses of the disk 
increase to O.IMq. 

Figure [T3] summarizes our results and shows the param- 
eters (orientation and magnetization) for which disks can 
form. 



6.3. Velocity 

To estimate the rotational velocity in the disk, we average 
v^{r, (J3, z) azimuthally and axially over the thickness of the 
disk. We can compare it to the Keplerian velocity VK(r) = 
y^GM{r)/r, with M{r) the mass within a sphere of radius 
r. We find that the rotational velocity in the disk is nearly 
Keplerian, as shown in Fig. [T5]for fi = 5, a — 45 and 90°, 
and /i = 3, a = 90°. In the perpendicular case, for both 
magnetizations, the rotation velocity is nonetheless slightly 
sub-Keplerian; for a = 45°, the disk has a fiat rotation 
curve. 

To understand the presence of sub-Keplerian rota- 
tion velocity profiles, we plot in Figure [11] the ratio of 
the magnetic pressure to the rotational kinetic energy 
\ogiQ{PB/Ek) — \ogiQ{B^/4:Trpv'^) for the same magnetiza- 
tions and angles. This shows that these velocities are sub- 
Keplerian because of the magnetic support in the disk: the 
more tilted the axis of rotation, the lower the ratio of the 
magnetic pressure to the rotational kinetic energy. In par- 
ticular, in the perpendicular case, the vertical component of 
the magnetic field is smaller, whereas the density is higher 
since the disk is more massive. 



6.4. Disk shape 

To estimate the radius and height of the disks, we use 
the first four criteria defined in section 16.11 (c/. Fig. [T^ . 
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Figu re 13. Mass of the di sk as a function of tim e for ji — 17 
(Fig. [13R1), fi ^ 5 (Fig. pM, ^ = 3 (Fig. mm, and 
fi = i (yig.|l3(d)|. 



However, as shown below, using only these criteria to in- 
fer the disk radius and height leads to an overestimate of 
their values, because of the large size of the spiral arms as 
evidently seen in Fig. 1171 We thus consider them as upper 
limits for the disk radius and height. For a more realistic 
estimate, we constrain the disk to the central, denser ro- 
tating object and therefore consider the region of the disk 
with a density above 10^ cm~'^, which is our last criterion. 
When the disk forms, its radius generally evolves with 
time, reaching values between 200 and 400 AU (for the min- 
imum estimate of the radius, -Rmin) from 500 to 800 AU (for 
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Figure 14. Disk formation in the parameter space investi- 
gated by the simulations (incUnation angle a versus magne- 
tization /i). The red diamonds are configurations in which 
approximately Keplerian disks form, the grey squares de- 
note configurations in which disks form with flat rotation 
curve (c/. Fig. [T5l and the related discussion for more de- 
tails); the white circles are configurations with no signifi- 
cant disk (Mdisk < S.fO-^Mg). 



the maximum estimate, i?max)- However, much smaller val- 
ues are obtained when the braking is strong. The minimum 
height i/min is about from 20 to 40 AU, while the maximum 
height -ffmax reaches 140 AU. These values can be compared 
to analytical estimates of the characteristic height of a hy- 
drostatic disk, i?th — \/cs/47rG'/9, taking for Cg and p their 
mean values in the disk. The lower height estimates are in 
good agreement with the theoretical values. All results are 
summarized in Table [51 Finally, Figure [T7] shows two den- 
sity slices (at the beginning and the end of the simulation) 
in the perpendicular case for ^ = 5. A central well-shaped 
disk and two large spiral arms are clearly visible. We note 
that in our estimates, the maximum estimated radius cor- 
responds to the maximum radius of the spiral arms. In con- 
trast, the minimum radius accuratly describes the central 
denser object. 

Figures [13 [TH and [ID] show, in contrast to Fig. [T71 
density slices in t he equ atori al pla ne when no massive 
disk forms. In Fig. 18(a) and 19(a) the pseudo-disk can 



be clearly distinguished around the protostellar core. In 
Fig. 1201 the two arms correspond to matter collapsing along 
the magnetic field lines. 

6.5. Discussions 

6.5.1. Impact of the criteria 

Estimates of the disk mass generally depend on the working 
definition of the disk, and may lead to large overestimates. 
One example is to calculate the disk mass using a simple 
criterion based on a comparison between rotation and infall 
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Figure 15. Radial profile of rotational velocity in the disk, 
for /i = 5, a = 45 and 90° and /i = 3, a = 90°. The straight 
lines are the rotational velocity for different time-steps; the 
dotted lines are the Keplerian velocity. 



ve locity 



Vd, > Vr, wh i ch is t he one used for example 



in iMachida fc Matsumotol (|201l[ )). The panels in Fig. [21] 
display the mass of the "disk" found with this criterion, 
which range from 0.3 to 0.5 Mq in all cases. They are more 
massive in the more tilted cases than in the aligned one, 
and for higher magnetizations, they are less massive, even if 
their formation is not prevented. Thus, while the trends are 
similar to those that we found previously, the mass of the 
disk can be greatly overestimated by using such a criterion. 

6.5.2. Comparison with observations 

Several studies have tried to infer disk masses from low 
resolution observations (with spatial resoluti ons of about 
250 AU), without resolving the disk itself (JEnoch et al.l 
I2OO9L [2OIII) ■ For this purpose, they used a detailed emission 
model, coupled with an analytical model for the envelope, 
the cavity, and the disk. The envelope model is that of a ro- 
tating, collapsing sphere developed bv lUlrichI (|l976f ). with 
a cavity in which the density is set to zero to mimic the 
outflows. The disk density is given by a power-law depen- 
dence in rad i us an d a Gaussian dependence in height (see 
[Enoch et al.l (|2009f ) for more details). Using these profiles, 
they ran a grid of models to find the parameters that most 
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Table 2. Disks characteristics (maximum mass of the protostellar core (which corresponds to M{n > 10 cm ^)) 
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Figure 16. Ratio of magnetic pressure to rotational kinetic 
energy in the disk, for /i = 5, a = 45 and 90° and /i = 
3,Q = 90°. 



closely fit their observations, by computing detailed radia- 



tive transfer. Their most important parameters are the disk 
mass and radius. 

Following the same idea, we attempt to deduce the disk 
mass from our simulations using their method, but without 
the radiative transfer. At several time-steps, we compute 
column-density maps of our simulations, with an angle be- 
tween the axis of rotation of the core and the line of sight 
of 15° (which corresponds to the angle of the line of sight 
in their best- fit model) . We consider the total mass of our 
simulation, IMq, to be the mass of the envelope, and adopt 
the same power-laws f or the density profiles of the env elope 
and the disk, namely (|Ulrichlll97a lEnoch et al.ll2009f ) 

(28) 

(29) 

^/"^ the height of the 



Pcnv 
Pdisk 

with 



cx r 



ex r 



-1.5 



-(r/H{r)f 



H{r) = r{Ho/Rdisk){r/Rdisi. 
disk. We take an outfiow opening angle of 20° (which is 
the one tha t most closely fits their observations) . As in 
lEnoch et al.l (|2009D . we infer a centrifugal radius Re from 
column-density profile corresponding to the radius where 
the slope of the column-density profile changes. The disk 
radius is another parameter that is varied; for simplic- 
ity, we only consider two different disk radii, i?disk = Re 
and Rc/'2. The disk vertical scale-height Hq is equal to 
0.2i?disk- With these parameters, we run the models with 
Mdisk = 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5 Mq. The model 
grid, which is a grid of column-density maps, is compared 
to snapshots of our simulations by means of a mean squared 
error analysis to find the best-fit model. 

Using this method, we find that the best fit to our sim- 
ulations is characterized by a disk radius -Rdisk = -Rc/2 and 
masses between 0.4 and 0.5 AIq for ^ = 5, masses between 
0.1 and 0.5 Mq for ^ — 3, and masses between 0.2 and 0.5 
Mq ioT h = 2 (see Fig. l^ . 

When massive disks actually form, which is not the case 
for the stronger magnetizations, this method provides re- 
sults that are in relatively good agreement with the disk 
masses inferred from our simulations (with an overestimate 
of 30 to 40 %). However, this comparison shows that in 
most cases density structures are detected that are mis- 
taken for disks, particularly in either the aligned case or 
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Figure 17. Density slice in the region of the disk for /i 
a = 90° at t = 20 000 and 25 000 yr. 
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Figure 18. Density sHce in the equatorial plane for fi — 5, 
a = 0° at t = 20 000 and 25 000 yr. 



for the higher magnetizations when no disks form, leading 
to a large overestimate of disk masses. Those density struc- 
tures, which also include the cavity of the outflows, actually 
correspond to a selection effect resulting from the use of a 
simple velocity criterion {v^ > Vr). 

To verify this assumption, we remove the cavity of the 
outflows in the /i = 5, a = 0° case and repeat the analy- 
sis. Figure [231 shows column-density maps after taking into 
account all the gas (left panel) and removing the gas that 
belongs to the cavity (right panel). In the central region, 
within 500 AU, there is a discrepancy in the column-density 
of a factor three between those two maps: the cavity is a 
massive structure, which can be mistaken for a disk by pro- 
jection effects. The best-fit model for a snapshot of this 
simulation without the cavity is Afdisk — O.Olil/0, where 
it was Mdisk = O.5M0 for the simulation with the cavity. 
Therefore, observed massive disks may actually be outflow 
cavities. 



6.5.3. Time formation of disks 

Our study shows that in very magnetized cases (^ ^ 1 — 3) 
disks beginning to form at the earliest time of star forma- 
tion (in Class stage) will not eventually form or remain 



small. This may appear to contradict the ubiquity of disk s 
at the Class I and later phases (e.g. iHaisch et al.l 12001"). 
However, it is worth stressing that the magnetic braking 
represents an exchange of angular momentum between the 
inner and outer parts of the prestellar core. In particular, 
the envelope should then operate as a reservoir able to ac- 
cept the excess of angular momentum present in the dens- 
est regions of the prestellar core. As accretion proceeds, the 
mass is the envelope diminishes and it is unlikely that mag- 
netic braking remains efficient. We therefore speculate that 
disks will always form but that their formation time and 
size will depend strongly on the magnetization and the an- 
gle between the rotation axis and the magnetic field; later 
on, we may be able to reduce the formation time and in- 
crease the size. 

Our conclusion - that both the formation time and disk 
size depend stro ngly on the magnetiza tion - is qualitatively 
similar to that of lDapp fc Basul (J2010r ) , although the under- 
lying reason is different. 



7. Conclusions 

We have presented an analytical analysis of a collapsing 
magnetic cloud that demonstrates that the magnetic field 
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We have then presented simulations of the collapse of 
can remove angular momentum less efficiently when the prestellar dense cores with different magnetizations /i and 
rotation axis is perpendicular to the magnetic field than in both aligned and various misaligned configurations. The 
when they are both aligned. orientation of the rotation axis with respect to the mag- 
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Figure 21. Mass o f the " disk" evolutio n for / i 
2TR1), A^^5 (Fig. [2T(b)] ), M 



17 (Fig. 
3 (Fig.[2T(c)|), and ^ = 2 



21(d)), obtained with a very simple rotation criteria 



(w0 > Vr). The masses are significantly overestimated. 



netic field, a, has a strong effect on the formation of the 
adiabatic first core and the disk formation. In particular, 
we have performed a detailed analysis of the transport of 
the angular momentum in the simulations, and character- 
ized the disks when they formed. Our main results are the 
following: 
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Figure 22. Obse rvation al estimate of t he mas s of the "disk" 
for ^ = 5 (Fig. [22(a)| ), ^ = 3 (Fig. [22(b)| ), and /i = 2 
(Fig. 22(c) ). These masses are inferred using a comparison 
between our simulations and analytical density profiles for 
the envelope, the outfiows, and the disk (see text). 



Misalignment has a strong impact on the outfiows and 
can suppress them; consequently, the angular momen- 
tum transport by the outflows decreases with a. 
Angular momentum transport by gravity increases with 
a, owing to the presence of the disks, particularly their 
asymmetric structures. 
The mass in the disks increases with a. 
For increasing magnetic fields, the disk masses decrease, 
with a limiting case being that of /i = 2, where disk 
formation is prevented. 

disks have typical mass up to 0.3 Mq and typical radii 
of from about 200 to 400 AU. 

In general, magnetic braking is the most important 
mechanism for transporting angular momentum. It al- 
ways dominates the transport of angular momentum by 
the fiow and, except for low magnetization (^ > 17), 
also dominates the transport by means of gravitational 
torques. 



— Magnetic braking decreases with a, but increases with 



We have shown that our conclusions depends on the crite- 
rion we use to define disks. We conclude that a simple ro- 
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tation criteria is insufficient and leads to estimates of disk 
masses that are far too high. 

We also analyzed our s i mulat ions following the method 
described in lEnoch et alJ (|2009f ). demonstrating that low 
resolution observations can mistake density structures for 
disks. 
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Appendix A: Euler's equation and angular 
momentum transports 

The Euler's equation for a magnetized fluid can be written 
as 

pa,v + p(v.V)v = -v(^P+|^)-pg+(^^.v)B, 

with p the density, v the velocity, P the gas pressure, B the 
magnetic field, and g the gravitational acceleration where 
g = V<i> and $ is the gravitational potential. 

In cylindrical coordinates, the azimuthal component of 
the left hand side of lA. II can be written 

pdtVij, + p {vrdrv^ + — 1 d^v^ + v^dzvA . (A. 2) 

Using the continuity equation 

V • (pv) +dtp ^ Q (A.3) 

-dr{rpvr) + -dtj,{pv^) -f dyXpVz) + dtp = 0, (A. 4) 

it can be written as 
pdtv^ -f v^dtp + p—dr{rv^) H dr{rpvr) 

+ p — ^^(w^)-! d^{pv^) 

+ pdzV^ +v^dz{pvz) 

= dtipv^) + -dr{rpVrvA + p— — 
r r 

1 



-d^{pv^v^) + dzipvzv^) 



dtipv^) + V • (pw^v) -f- p 



VrV^ 



(A.5) 



We can do the same for the magnetic tension component 
of the equation 

^ (srdrB^ + ^^ + ^d^B^ + BAB^ , (A.6) 

using the solenoidal constraint (V • B = 0); it comes 

_L(v.(B,B, + Mi). ,A.7) 



The density p can be expressed as a function of the 
gravitational acceleration, using the Poisson's equation 



1 



AttG 



Vg. 



(A.8) 



The azimuthal component of the gravitational term of 
Euler's equation thus becomes 



P50 



50 



Vg 



47rG 
AttG \r 



dr{rgr) + -d^g^ + dz 

r 



1 



AttG 



— 1- dr{grg<t,) + -905030 + 9zgz 

r 1- J. f 



1 



grdrg4> + -g^d^g^ + gzdzg^, 



(A.9) 



and using Schwarz' theorem 
1 



P94> 



pg<t> 



AttG 



'-^^ + -dr{rgrg4>) + -d^g^g^ + dzgzg<i> 



^(^ + V.(50g)-i(e0.V)g^), (A.IO) 



we can identify a curvature term {grg4>/f)i a tension term 
(V • (ff0g)), and a pressure term ((e0 • V) g^/8nG). 

Therefore, the azimuthal component of the Euler's 
equation can eventually be written as 



dt{pv^) + V • (/9U0v) + p 



VrV^ 



B^ .g2 



-(^^•^)l^' 8. SnG 



+ ^ ( V . (i30B) ' ^^ 



47r 
1 



^g) + 



r 
9rg4, 



inG V '"' "' r 

If we multiply this equation by r, it becomes 

B^_J^\ 



(A.ll) 



dt (prv^) 



V -r 


pv^v + 


^^B 

An 


+ ^" gl 



^+8;r-8^r^ 



0, 



(A.12) 



since r is time-invariant. This equation expresses the an- 
gular momentum conservation; we can identify the mag- 
netic and the gravitational torques V • (ri?0B/47r) and 
V • (rgcjig/ATTG) with i G {r,z}, which are responsible for 
the angular momentum transport by means of a magnetic 
field and gravitation, respectively. 



Appendix B: Convergence 

Additional sets of simulations were run to test the numer- 
ical convergence. These sets of simulations show that the 
numerical dissipation does not significantly change our re- 
sults, which are qualitatively invariant. 

In the first set of simulations, we change the Jeans re- 
finement strategy to increase the spatial resolution. A cell 
was previously refined if its size exceeded one-tenth of a 
Jeans' length (csin/Gp)^'"^). We run simulations with 15 
cells per Jeans' length (HRl) and another one with 20 cells 
per Jeans' length (HR2). 
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In the last set of simulations, we change the Courant 
number (from 0.8 to 0.4) to increase the temporal resolu- 
tion. 

All sets were run with /x = 5 and 2, and a — 0, 45, and 
90° for increasing spatial resolution, and with ^ — 5, and 
a = 0,45, and 90° for increasing temporal resolution. 

Our results remain qualitatively similar even if conver- 
gence is not achieved. These convergence runs show that the 
mass of the disk is slightly overestimated in our previous 
analysis. 



B.l. High spatial resolution simulations 



U=5 



B.1.1. Angular momentum 



Figures IB. II show the evolution of specific angular momen- 
tum for n > 10^ cm~^ for, respectively, fj, — 5 and fi — 2. 
The left panel display its evolution for 10 points, the cen- 
tral panel for 15 points, and the right panel for 20 points 
resolved per Jeans' length. For all angles {a = 0, 45 and 
90°); the specific angular momentum decreases with the 
resolution. 

With increasing resolution, more momentum is trans- 
ported by the magnetic field and outflows, and the steady 
state is reached at a later stage. 

It is clear that numerical convergence has not yet 
been reached and that treating magnetic braking requires 
a very high spatial r esolution (iCommerqon et all l2010t 
iHennebelle et al.ll2011ll . 

Even though these results are quantitatively different 
from the simulations presented in this paper, they are qual- 
itatively similar: magnetic braking still decreases signifi- 
cantly when the angle between the rotation axis and the 
magnetic field increases, and a comparable amount of an- 
gular momentum is carried away by magnetic braking and 
outflows. 



B.l. 2. Disk mass 

The mass of the disk stays roughly the same or decreases 
with increasing resolution, as shown in Fig. IB. 21 The mass 
of the disk can decrease with increasing resolution because 
of the more efficient transport of angular momentum in the 
higher resolution cases. For example for ^ = 5, a = 90°, the 
mass of the disk at t = 25 500 yr is 0.18 Mq (LR) compared 
to 0.15 Mq (HRl) and 0.13 Mq (HR2). 

However, our previous conclusions still hold: when mag- 
netization is relatively strong (/x = 5), disks form only when 
the rotation axis is misaligned with the magnetic field, and 
for lower /x (meaning stronger magnetization), magnetic 
braking acts so strongly it prevents disk formation. 



B.2. High temporal resolution simulations 
B.2.1. Angular momentum 



Figure IB. 31 displays the evolution of the specific angular 
momentum for n > 10^ cin~'^, /i = 5, for a Courant number 
of 0.8 (left panel) and 0.4 (right panel), which corresponds 
to a smaller time-step. As the resolution increases, more 
momentum is transported outward. 




19000 20000 



21000 22000 
time (yr) 



23000 24000 25000 



Figure B. 4. Disk mass for /x = 5 and a Courant number 
of 0.4. 



B.2.2. Disk mass 



Figure IB. 41 shows the evolution of the mass of the disk for 
Courant number of 0.4. 
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Figure B.l. Specific angular momentnm for /j, = 5 (upper panels) and /i = 2(lower panels), n > 10^ cm '^, for a ~ 0,45, 
and 90°, with 10, 15, and 20 resolved points per Jeans' length. 
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Figure B.2. Mass of the disk as a function of time for /^ = 5 (upper panels) and fi — 2 (lower panels), for a = 0, 45, and 
90°, with 10, 15, and 20 resolved points per Jeans' length. 
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Figure B. 3. Specific angular momentum ior fi — 5,n > 10^ cm ■^j with a Courant number of 0.8 and 0.4. 
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